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Abstract. We calculate the meson screening mass in a pion superfluid in the framework of Nambu-Jona- 
Lasinio model. The minimum of the attractive quark potential is always located at the phase boundary of 
pion superfluid. Different from the temperature and baryon density effect, the potential at finite isospin 
density can not be efficiently suppressed and the matter is always in a strongly coupled phase due to the 
Goldstone mode in the pion superfluid. 



1 Introduction 

Understanding the behavior of quantum chromodynam- 
ics (QCD) at finite temperature and density is essential 
for a description of the development of the early universe 
and compact stars and for the explanation of the results 
from high energy nuclear collisions at Relativistic Heavy 
Ion Collider (RHIC) and Large Hadron Collider (LHC), 
where high temperature and/or density can be reached. 
From Lattice simulations QUEEE at finite temperature, 
there are two QCD phase transitions, namely the decon- 
finement transition from hadronic matter to quark matter 
and the chiral transition from chiral symmetry breaking 
to its restoration. 

At finite baryon density, the QCD phase structure be- 
comes more rich. Besides the deconfinement and chiral 
phase transitions, the color symmetry which is strict in 
vacuum is spontaneously broken and leads to a new phase 
at high density, the so-called color superconductor [S1[SJ[7] 
which may exist in the core of dense quark stars. How- 
ever, due to the fermion sign problem [8], there is not yet 
precise lattice result at finite baryon density. While the 
properties of the new phases at extremely high density 
can be studied with perturbativc QCD, to investigate the 
phase transitions themselves at moderate baryon density 
one needs effective QCD models at low energy. 

The isospin symmetry is spontaneously broken when 
the isospin chemical potential is larger than the pion mass, 
resulting in the so-called pion superfluid phase [9] which 
may exist in the core of neutron stars. In comparison with 
the temperature and baryon density effect, the study on 
isospin effect or pion superfluid phase includes new phe- 
nomena: 1) The isospin structure of QCD phase transi- 
tions can be directly examined through lattice simulations 
without serious technical problems [TOllTTj emerged at fi- 
nite baryon density; 2) While for the color superconduc- 
tivity at finite baryon density the Goldstone modes corre- 
sponding to the local color symmetry breaking are eaten 
up by the Higgs mechanism, the Goldstone mode corre- 
sponding to the global isospin symmetry breaking dom- 
inates the thermal and dynamic properties of the pion 



superfluid [T2J; 3) While the phase transitions of decon- 
finement and chiral symmetry restoration occur at finite 
temperature and baryon density, there might be no decon- 
finement |13j and the chiral structure is also significantly 
changed [14] at finite isospin density. 

From the lattice calculated thermodynamics at finite 
temperature, the energy density [3] and pressure U do 
not reach the corresponding ideal gas limit in the decon- 
fincd phase, this means that the quark matter close to the 
phase transition is not a weakly but a strongly interact- 
ing system. The strongly coupled quark matter predicted 
by the theoretical calculations is recently supported by 
the strong collective flow [tTniTMTT] observed in heavy ion 
collisions at RHIC. It is widely accepted that the coupling 
strength will decrease with increasing temperature and 
baryon density, and the quark matter will be in a weakly 
coupled phase when the temperature or density is high 
enough. 

The potential between two quarks in hot and dense 
medium can well describe the changes in quark proper- 
ties at finite temperature and density. It presents a di- 
rect way to understand whether the matter is a strongly 
coupled one. For instance, the heavy quark potential ex- 
tracted from lattice simulations [IB] at finite temperature 
is the key point to explain the J/tp suppression [T2] ob- 
served in heavy ion collisions at RHIC. At finite temper- 
ature and baryon density, the effect of chiral symmetry 
restoration on the quark potential is investigated |20) in 
the framework of Nambu-Jona-Lasinio (NJL) model |21) . 
It is found that the minimum of the attractive potential 
is located at the critical temperature T c or the baryon 
chemical potential fi c B of the chiral phase transition, and 
the strongly coupled matter survives only in the temper- 
ature region 1 < T/T c < 2 — 3 or baryon density region 
1 < hb/h% < 2. In this paper, we study the meson screen- 
ing mass and the quark potential at finite isospin density 
in the NJL model. We will focus on the potential in the 
pion superfluid and explain its surprising isospin behav- 
ior induced by the Goldstone mode corresponding to the 
spontaneous isospin symmetry breaking. 
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The original NJL model [2T] is inspired by the BCS 
theory describing normal electron superconductor, and there- 
fore its version at quark level is widely used to describe chi- 
ral symmetry restoration, color superconductor and pion 
superfluid (For reviews and general references, see Ref. [2"2l 
To understand the results of lattice 
QCD thermodynamics in terms of quasi-particle degrees of 
freedom, the model is recently extended to include Polyakov 
loop dynamics (PNJL) - Because of the contact 

interaction among quarks, there is no confinement in the 
model and it is necessary to introduce a momentum cut- 
off A to avoid the singularity in momentum integrations. 
From the uncertainty principle, the minimal length scale 
in the NJL model is R ~ 1/A. Taking the standard cut- 
off value A ~ 600 MeV \gZ$EMMKMM > we have R ~ V 3 
fm. This means that, the quark potential calculated in the 
NJL model might be reasonable in the region of r > 1/3 
fm. In the short range the model is probably not applica- 
ble. 

The paper is organized as follows. In Section [5] we re- 
view the meson propagators in the pion superfluid in the 
NJL model and calculate the meson screening masses. In 
Section [3] we calculate the quark potential exactly and 
in pole approximation, analyze the dominant contribution 
from the Goldstone mode due to spontaneous isospin sym- 
metry breaking, and discuss the strongly coupled matter 
in the whole pion superfluid. We summarize in Section [H 



2 meson screening masses 

To simplify the numerical calculations, we take the orig- 
inal NJL model and neglect the Polyakov loop poten- 
tial [29 , 30 , 31 which disappears at zero temperature. The 
two flavor NJL model is defined through the Lagrangian 
density 

r 2 2" 

C = ip ii^d u - m + fijo) ip + G (ipip) + (t/ji^ri/j) 

with scalar and pscudoscalar interactions corresponding to 
a and n excitations, where mo is the current quark mass, 
G is the coupling constant with dimension (GcV)~ 2 , the 
quark fields ip and ?/>, the Pauli operators r = (ti,T2,ts) 
and the quark chemical potential fi = diag{fi u , fid) = 
diag(fiB/3 + \ii /2, fis/3 — P-i/2) are matrices defined in 
flavor space, and fiB and \i\ are baryon and isospin chem- 
ical potentials. At /i/ = 0, the system has the symme- 
try J7b(1) <S> SUi{2) (g) SUa(2), corresponding to baryon 
number symmetry, isospin symmetry and chiral symme- 
try. However, the isospin symmetry SUj(2) is explicitly 
broken down to Uj(l) global symmetry at small /ij, and 
then the Ui(l) symmetry is further spontaneously broken 
down with condensation of charged pions at large p,i . At 
[lb = 0, the Fermi surfaces of u(d) and anti-d(u) quarks 
coincide and hence the condensate of u and anti-d quarks 
is favored at /ij > and the condensate of d and anti-M 
quarks is favored at fii < 0. A finite [Ib provides a mis- 
match between the two Fermi surfaces and will reduce the 
pion condensation. 

Introducing the chiral condensate 

a = {HA (2) 



and pion condensate 

7r = v / 2(-0i75T + V) = V / 2(V'i75T_V') (3) 

with t± = (ti ± it-i) /\2, which are respectively the order 
parameters for the chiral phase transition and pion super- 
fluid, the inverse of the quark propagator in mean field 
approximation defined in flavor space 



Sip) 



Suuip) Sudijp) 

Sdu{p) Sdd(p) 



can be written as 



2iGn l5 



m 2iGir~f5 
7 M £V + Vdjo - m 



(4) 



(5) 



where m = mo — 2Ga is the dynamic quark mass induced 
by the spontaneous chiral symmetry breaking. With the 
quark propagator, the physical condensates a and tt are 
determined by the gap equations, 



d 4 p 



(2nY 



Tr [S uu (j>) +Sddip)} , 



d 4 p 



Tr [iSudip) + Sdu{p)) 75] 



(0) 



with the trace taken in color and Dirac spaces, where 
the four dimensional momentum integration is defined as 
/ (i 4 p/(27r) 4 = iT J d 3 p/{2n) 3 J2 n ^ n * ne imaginary time 
formalism of finite temperature field theory with the Mat- 
subara frequency w„ = (2n + VjitnT, n = 0, ±1, ±2, ... for 
fermions. 

In the NJL model, the meson modes are regarded as 
quantum fluctuations above the mean field and can be ef- 
fectively expressed at quark level in terms of quark bubble 
summation in random phase approximation (RPA) |22[ 
23ll241l25ll26] . In normal phase without pion condensation, 
the bubble summation selects its specific isospin channel 
by choosing at each stage the same proper polarization. 
In the pion superfluid phase, however, the quark prop- 
agator contains off-diagonal elements in flavor space, we 
must consider all possible isospin channels in the bubble 
summation. In this case, the meson polarization function 
becomes a matrix in the four dimensional meson space 
with off diagonal elements [TJ] 



n 3 kig) 



d A p 

J2n) 4 



Tr 



r*s 



p 



as P 



(J) 



with j, k = a, 7r + , 7r_, 7To, where the trace is taken in 
color, flavor and Dirac spaces, and the meson vertexes are 
defined as 



r, 



1 3 


= a 


iT +l5 j 


= 7T+ 


iT-75 j 


= 7T_ 


ir 3 j5 j 


= 7T 



1 j =(T 

*T-75 j = 7T+ 
iT +l5 j = 7T- 
«T"375 j = 7T . 



(8) 



From the definition (|7|) , the polarization function II is 
a symmetric matrix in the meson space with II jk = H kj , 
and the neutral pion 7Tq is decoupled from the charged 
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pions 7r + and 7r_ and the isospin singlet cr with 77„- 0(T = 
n iT0 T T+ = iI 7ro7r _ = 0. Considering further the relations 



77, 



T _{q) = n a7r+ (-q) and II n _ n _(q) = II v+v+ (-q), there 



are only five independent polarization elements n aa , IJ, 



77 lr+lr+ , II 7r+7l _ and II no7ro . Taking the trace and perform- 
ing the fermion frequency summation in ([7]), and intro- 
ducing the following definitions based on the quark energy 



(E p + 6/X//2) 2 + 4G 2 tt 2 



,ab 



ab 



v /(l + (^ + M //2) /ef) /2, E a p = aE p , vf - 
with a, b = ± and the Fermi-Dirac distribution function 
/(a;) = 1/ (e^ x ~ IJ ' B ' s ^' T + l), the complicated expressions 
for the polarizations are greatly simplified, and the five in- 
dependent elements can be explicitly expressed in a com- 
pact way in (9), with p± = p ± q/2. Note that there is 
no Lorentz invariance at finite temperature and density, 
and the polarization elements are no longer functions of 
the Lorentz scalar quantity q 2 = q 2 — q 2 but depend sep- 
arately on <7o and q 2 . 

Describing the meson exchange in quark scattering via 
the quark bubble summation and taking into account the 
spontaneous isospin symmetry breaking, the meson prop- 
agator can effectively be expressed as a matrix in the four 
dimensional meson space |14) 



£(<?o,q 2 



2G 



1-2G77 (^q 2 )' 



(10) 



The dynamical meson mass Mj in vacuum is defined as the 
pole of the meson propagator at q 2 = M 2 . At finite tem- 
perature and density, the position of the pole is changed 
A7 2 ,q 2 =0, 



to q 2 



det[l-2G77(M 2 ,0)] = 0, 



(11) 



where the symbol det means the determinant of the ma- 
trix 1 — 2G77. In normal phase without pion condensa- 
tion, it is simplified to four independent mass equations, 



1 - 2G7L- 



0) = 0. At fxi = and in chiral limit, 



it is found [32] that the three pions are the Goldstone 
modes corresponding to the spontaneous chiral symmetry 



det 



1 



breaking and dominate the thermodynamics of the sys- 
tem at low T and low [Ib- In the pion superfluid, how- 
ever, mesons a, n+ and 7r_ are coupled together and no 
longer the eigen modes of the Hamiltonian of the system, 
the new eigen modes, labeled by cf, 7f+ and W _ are linear 
combinations of cr, 7r + and 7r_ . The pole equation in 
this case is separated into 1 — 2G77 7ro7ro (M| , 0) = for 
7Fo (to avoid confusion with the itq in normal phase, we 
still use 7fo to stand for ttq in the pion superfluid) and 

2G77 (M 2 , 0) = for W,W+ and tt_, where the 

polarization 77 is the sub-matrix of 77 defined in the me- 
son subspacc {cr, 7r + , 7r_}. The T, (1b and /i/ dependence 
of Mj is shown in Ref . [33] . 

From the Yukawa potential between two nucleons, V(r) 
e~ Mr /r, its strength is governed by the mass M. of the 
exchanged boson. The long range force is suppressed or 
screened by the massive boson (this is the reason why M. 
is called the screening mass). Its inverse is the screening 
length 7Z = 1/Ai. At finite temperature, the screening 
mass is defined as the pole of the boson propagator at 
g 2 = 0, q 2 = — M 2 . In the NJL model at quark level, the 
meson screening mass A4 which control the quark poten- 
tial arc defined as 



det [l-2G77(0,-A4 2 )] = 0, 



(12) 



which is again simplified to four independent equations 1 — 
2G77jj(0, — -M 2 ) = in normal matter for j = <t,tt + ,tt^,ttq 
and separated into 1 — 2GII 7To7To (0, — ) = for tt and 

det[l — 2G77(0, — M 2 )] — in the meson subspacc for 
cf, 7f+ and 7f_ in the pion superfluid. 

Before the numerical calculation, we can analytically 
prove that in the pion superfluid the mass equation (fTT]) 
at M = or (fl"2"|) at A4 = becomes exactly the second 
gap equation of (JB]) for the pion condensate, 



det [1 - 2G77 (0, 0)] = 0. 



(13) 



This means that the Goldstone mode with M = M. = 
0, corresponding to the spontaneous isospin symmetry 
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breaking, is automatically guaranteed in the RPA approx- 
imation in the NJL model, and the Yukawa potential via 
the exchange of the Goldstone meson is not screened. 

There are three parameters in the NJL model, the cur- 
rent quark mass mo, the coupling constant G and the 
momentum cutoff A. In the following numerical calcula- 
tions, we take m =5 MeV, G=4.93 GcV -2 and yl=653 
McV. This group of parameters corresponds to the pion 
mass M !r =134 MeV, the pion decay constant .^=93 MeV 
and the effective quark mass m=310 MeV in vacuum. 
To extract reliable conclusions from the numerical cal- 
culations in the model, the thermodynamic parameters, 
namely the temperature and quark chemical potentials, 
should be much less than the cutoff, T, jUg/3, /if/2 <c A. 
In our following numerical calculations we take T < 300 
MeV, he < 900 MeV and in < 600 MeV. Considering 
the critical isospin chemical potential fij — 134 MeV [T4"] . 
\ii = 600 MeV ~ 4fij is large enough for the discussion of 
pion superfluid. 

With the coupled gap equations © for the chiral and 
pion condensates, the phase boundary of the pion super- 
fluid as a hypersurface in the three dimensional space of 
temperature and baryon and isospin chemical potentials is 
shown in Fig. [1] Since the pion condensate increases with 
decreasing temperature and baryon density and increas- 
ing isospin density, the pion superfluid phase is located 
in the region of low temperature, low baryon density and 
high isospin density. By comparing the second gap equa- 
tion of © for the pion condensate with the dynamical 
mass equation (TTTj) for pions, the critical isospin chemical 
potential for the pion superfluid at T = Hb = is exactly 
the pion mass in vacuum, /ij = M n = 134 McV. This 
critical value increases with increasing temperature and 
baryon chemical potential. 




Fig. 1. The phase transition hypersurface of pion superfluid 
in the three dimensional space of temperature T and baryon 
and isospin chemical potentials hb and fii in the NJL model. 

The isospin chemical potential dependence of the screen- 
ing masses at fixed temperature and baryon chemical po- 
tential is shown in Fig. f2J For the chosen T = 100 McV 



I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 




0.0 0.1 0.2 0.3 0.4 

H, (GeV) 

Fig. 2. The meson screening masses M as functions of isospin 
chemical potential at temperature T = 100 MeV and baryon 
chemical potential — 300 MeV. 



and fJ-B = 300 MeV, the critical \xi for the pion super- 
fluid is fij = 142 MeV. Since mesons carry isospin charge, 
the isospin dependence of the meson propagator D comes 
from not only the chiral and pion condensates cr(fii) and 
7r(/ij) but also an explicit shift for the meson energy qo — > 
q - M/, D((q - /4) 2 , q 2 ; d(/i/), tt(/i/)), where \i\ is the 
meson isospin chemical potential. This explicit \ii depen- 
dence leads to a dynamical mass splitting for the charged 
pions [14] even in the normal phase. For the screening 
mass defined at q 2 = 0, q 2 = — M 2 of the propagator, 
however, the explicit and inexplicit isospin dependence in 
the normal phase is the same for the two charged pions, 
and therefore they are degenerate in the normal phase. 
Considering the condensation of 7r + at > fij, 
and Aiw- are different from each other in the pion super- 
fluid. The degeneracy of M.-k + and M. n _ in the normal 
phase leads to two zero modes at the critical point, which 
is the reason why the strength of the quark potential ap- 
proaches to the maximum at the phase transition very 
fast, see the discussion below. In the pion superfluid, only 
ttq is still the eigen mode, but a, 7r + and 7r_ arc replaced 
by the new eigen modes a, 7r+ and 7F_ . 7f + is the Goldstone 
mode with zero screening mass Ai- !T+ =0, as we analyzed 
above. At small isospin chemical potential /i/ < /ij, the 
chiral symmetry S'L t a(2) is already explicitly broken to 
Ua(X)- Therefore, the chiral symmetry restoration at high 
isospin chemical potential means only degeneracy of a and 
7To, the charged 7r + and 7r_ behave differently. However, in 
Fig. [2] the degeneracy at high isospin chemical potential 
is not for 7T0 and W but for Wq and W _ . This is due to the 
strong mixing among a, 7r + and 7r_ |14| . 

In the top panel of Fig. [3] we show the screening masses 
in the normal phase as functions of temperature at fixed 
baryon and isospin chemical potentials [Ib = 300 MeV 
and [ii = 100 MeV. The charged pions 7r + and 7r_ coin- 
cide with each other in the whole plane, as we discussed 
above. Since there is no a — n mixing in the normal phase, 
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Fig. 3. The meson screening masses Ai as functions of tem- 
perature T in normal phase (top panel) at isospin chemical po- 
tential [ii — 100 MeV and in pion superfluid (bottom panel) at 
Hi = 200 MeV. In both phases the baryon chemical potential 
is fixed to be /i B = 300 MeV. 



the a and ttq mesons become degenerate at high temper- 
ature where the chiral symmetry is almost restored. For 
a small isospin chemical potential, the mass splitting be- 
tween neutral and charged pions is slight. 

In order to clearly see the effect of spontaneous and ex- 
plicit isospin symmetry breaking on the screening masses, 
we show in the bottom panel of Fig. [3] their temperature 
dependence at a higher isospin chemical potential. For the 
chosen \n = 200 MeV and [is = 300 MeV, the critical 
temperature for the phase transition from pion superfluid 
to normal phase is T c ~ 155 MeV. In the pion superfluid at 
T < T c , 7r + is the Goldstone boson with screening mass 
A4jf + = 0. Because of the mixing among er, n + and 7r_, 
7fo an d 7T-, instead of er and 7To, approach to each other 
at low temperature where the mixing is strong enough. 
When T is close to the critical point, W- approaches to 
the Goldstone mode due to the weakening of the pion 



condensation. In normal phase at T > T c , 7r + and 7r_ co- 
incide again, and with increasing temperature a and ttq 
mesons approach to each other when the chiral symmetry 
is gradually restored. 



3 quark potential 

By analogy with the Yukawa potential between two nu- 
cleons through one boson exchange, the static quark po- 
tential via one meson exchange can be expressed as the 
Fourier transform of the meson propagator at go = 0, 



V(r) 



■X>(0,q 2 



(14) 



In the flavor SU(2) NJL model, the meson propagator in 
the pion superfluid is a 4 x 4 matrix with off diagonal 
elements, the quark potential becomes 



V(r) 



d 3 q 

(2tt) 3 
G 



TrD(0,q 2 )e 4qr 
q s'm(qr) 



/ rf9Tr i-2Gi7(0,<? 2 )' 



(15) 



where the trace is taken in the meson space {er, 7r + , 7r_, ttq}. 
Note that there is no singularity for the propagator on the 
real momentum axis and the above integration is numer- 
ically straightforward without any technical problem. 

To better understand the potential through the screen- 
ing of the interaction, we take also the pole approximation 
in calculating the quark potential. In this case the poten- 
tial can be simplified as a summation over the poles, 



V(r) = - 



(27t) 3 



E 



■M) 



Airr 



-M, 



(16) 



with j = a - , 7r + , 7r_ , 7To in the normal phase and j = er, 7r + , H 
in the pion superfluid, where g.j q q is the effective meson- 
quark- antiquark coupling constant, 



2 

9jqq 



dot 1 



-2Gn(0,q 2 )] 



l-2G/7(0, 9 2 ) 



dq- 



■det [1 - 2G27(0,g 2 )] 



(17) 



-—Mi 



which is reduced to 



ijqq 



dn j3 (Q,q 2 



dq 2 



(18) 



q 2 = ~M2 



when the pion condensation disappears. 

In the normal phase at zero isospin density, it is easy to 
expect that the quark potential is gradually weaken by fi- 
nite temperature and baryon density, and this is confirmed 
by the lattice simulation [TH] and model calculation |20) . 
What we focus on in this paper is the isospin effect on the 
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Fig. 4. The quark potential as a function of the distance r 
between the two quarks in normal phase at hi — 100 MeV 
and in pion superfluid at fii = 200 MeV. In both phases the 
temperature and baryon chemical potential are fixed to be T — 
100 MeV and fi B = 300 MeV. The dotted line is the potential 
in vacuum. 



quark potential. In Fig. [4] we show the potential as a func- 
tion of the distance r between the two quarks in normal 
phase with hi = 100 MeV (dashed line) and in pion su- 
perfluid with iij = 200 MeV (solid line), the temperature 
and baryon chemical potential are fixed to be T = 100 
MeV and \ib = 300 MeV in the two phases. Consider- 
ing the momentum cutoff in the model, we do not discuss 
the potential in short distance and take the starting point 
r = 0.4 fm. The effect of the finite cutoff on the thermo- 
dynamics of the system is mainly at high temperature |32| . 
While this effect is not remarkable in the temperature and 
density region wc considered here, the potential is normal- 
ized to limr-j-oo y( r ) — in the high temperature limit. It 
is easy to see that this normalization is automatically as- 
sumed in the pole approximation (fT6|) . In both the normal 
phase and pion superfluid, the potentials oscillate slightly 
at large distance. This is from the so called Friedel oscilla- 
tion [33] induced by the sharp Fermi surface at zero tem- 
perature which has been widely discussed in different mat- 
ters 34,35,36 . Since the Fermi surface is smeared at finite 
temperature, the oscillation is suppressed by temperature 
and will disappear when the temperature is high enough. 
The difference between the normal matter and pion su- 
perfluid can be understood clearly in the pole approxima- 
tion (|16p . From the screening masses shown in Figs. [5] and 
[3l a (<r) is always much heavier than the other mesons 
in the region we discussed, its contribution to the poten- 
tial is almost screened, and in the intermediate and large 
distance the potential is mainly from the pion exchange. 
From Fig|3J the pion screened masses at hi = 100 and 200 
MeV satisfy the relations M Vo < M.^ and M.^_ < ■ 
If we neglect 7r + and W+, the potential in the pion su- 
perfluid should be above the potential in normal matter. 
However, when 7T+ and 7f+ are included, the mass relation 
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Fig. 5. The quark potential calculated exactly and in pole 
approximation in normal matter (top panel) with hi = 100 
MeV and in pion superfluid (bottom panel) with hi = 200 
MeV. In both phases the temperature and baryon chemical 
potential are fixed to be T = 100 MeV and /i B = 300 MeV. 



dependence, and the attractive potential in the pion su- 
perfluid may be stronger than the one in normal matter, 
namely the potential with hi = 200 MeV may be below 
the potential with hi = 100 MeV. This is confirmed from 
our numerical calculation, as shown in Fig. 2] While for 
Hi = 100 and 200 MeV the difference between the two 
potentials is small, the trend of the isospin dependence 
is really controlled by the Goldstone mode. We plotted 
also the potential in vacuum (dotted line) in FigJH To- 
tally different from the temperature and baryon density 
effect which weakens the interaction among quarks, the 
potential strength increases with isospin chemical poten- 
tial. 

The above numerical calculation for the potential is 
full but the analysis on its isospin dependence is based 
on the behavior of the screening masses. To know to what 
extent this analysis is valid, we compare in Fig.[5]the exact 
calculation with the pole approximation in normal matter 
with fii = 100 MeV and in pion superfluid with fii = 
200 MeV. Again the temperature and baryon chemical 
potential are fixed to be T = 100 MeV and hb — 300 
MeV in the two phases. While the pole approximation 
deviates from the full result remarkably in the region of 
r < 2 fm, and the Friedel oscillation is washed away in 
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the approximation, the screening masses do qualitatively 
describe the potential, for instance, the exponential decay 
and the saturation at large distance. 

To see clearly and comprehensively the medium effect, 
we show in Fig. [5] the isospin dependence of the quark 
potential at a fixed distance. We take r = 2 fm which is the 
typical diameter of a hadron, the potential at this distance 
can be used to describe the average medium effect on the 
coupling strength of the quark matter. From the phase 
diagram Fig. [TJ for the chosen baryon chemical potential 
\is = 300 MeV and the isospin region fij < 600 MeV, the 
maximum critical temperature for the pion superfluid is 
T c ~ 180 MeV. 




(GeV) 



Fig. 6. The quark potential at fixed distance r = 2 fm as a 
function of isospin chemical potential /ij at fixed baryon chem- 
ical potential fiB = 300 MeV and for different temperature T. 



Let us hrst look at a potential at a given temperature. 
For any T < T c , the strength of the attractive potential in- 
creases with increasing fij in normal phase and approaches 
to the maximum at the critical value fij. When the sys- 
tem enters the pion superfluid, the behavior of the poten- 
tial suddenly changes, its strength turns from increasing 
to decreasing and gets saturated fast. The sharp and deep 
valley at the critical point //j arises from the two zero 
modes 7r + and 7r_ or 7r + and 7f_ at the phase transition, 
shown in Figs. [2] and [3] Note that the saturated poten- 
tial in the pion superfluid is always below the vacuum 
potential V(r = 2 fm) ~ —0.06 GeV shown as dotted line 
in FigjSJ Therefore, from the valley structure of the po- 
tential, the quark matter is most strongly coupled at the 
phase transition point, and from the strong enough poten- 
tial at high isospin density, the pion superfluid is always 
strongly coupled. At T > T c , for instance T = 200 MeV, 
there is no more phase transition, and the valley structure 
disappears. At such high temperature, the potential at the 
distance of r = 2 fm almost disappears. 

We now focus on the temperature dependence of the 
potential at a given isospin chemical potential, shown in 




vacuum 

u= 100 MeV 

n 

H= 200 MeV 

|i=400 MeV 
|i=600 MeV 



0.0 0.1 0.2 

T (GeV) 



0.3 



Fig. 7. The quark potential at fixed distance r = 2 fm as a 
function of temperature T at fixed baryon chemical potential 
l_i B — 300 MeV and for different isospin chemical potential hi. 

Fig|7] For a small fij < 134 MeV, there is no pion su- 
perfluid even at zero temperature, the system is always in 
the normal phase. In this case, the temperature effect is 
under our expectation, the potential becomes weaker and 
weaker in the hot medium and finally vanishes. For a large 
Hi > 134 MeV, the system is in pion superfluid at T < T c 
and normal phase at T > T c . At a fixed isospin chemical 
potential, there is also a valley structure around the crit- 
ical temperature: The potential strength increases with 
increasing temperature and reaches the maximum at T c , 
and after crossing the phase transition point, the poten- 
tial changes suddenly and turns to become more and more 
weak and finally disappears. Different from the isospin 
structure shown in Fig(SJ the potential goes to zero very 
fast in the normal phase. 



4 Conclusion 

We have investigated the meson screening masses and 
quark potential in a pion superfluid in the frame of NJL 
model. The Goldstone mode corresponding to the sponta- 
neous isospin symmetry breaking plays a significant role in 
the thermodynamics of the system. The minimum of the 
attractive quark potential is always located on the phase 
transition hypersurface in the three dimensional space of 
temperature and baryon and isospin chemical potentials. 
While at extremely high temperature and baryon density, 
quark matter is expected to be in a weakly coupled state, 
the pion superfluid is in a strongly coupled phase even the 
isospin density is extremely high. These surprising prop- 
erties of pion superfluid may help us to understand the 
isospin asymmetric matter, like the core of compact stars. 

In this paper, we considered quarks in mean field ap- 
proximation and mesons in RPA. When we go beyond the 
mean field for the pion condensate by taking into account 
the feedback from the mesons, the phase boundary will be 
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shifted and the quark potential in the region of low tem- 
perature and low densities will be modified, like the case 
of chiral phase transition [37] . However, the minimum of 
the potential at the phase transition and the non-zero po- 
tential at extremely high isospin density will remain, since 
they are controlled by the Goldstone mode. 
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